Limit cycle induced frequency locking in self-sustained current oscillations of 

superlattices 
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The ac response of self- sustained current oscillations (SSCOs) in GaAs/AlAs superlattices (SLs) 
is derived based on the deformation of a limit cycle under an external ac driving force. Frequency 
locking into an integer fraction of the ac frequency is obtained in a periodic response in which a 
limit cycle deforms either with or without a topological change. This frequency locking is robust 
against the ac bias because a limit cycle can adjust itself. The results are verified both numerically 
and experimentally, indicating that SSCOs in SLs can be understood within the framework of the 
general concepts and principles of nonlinear physics. 
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Following the early pioneering studyP, 01 on vertical 
electron transport in superlattices (SLs), one of the re- 
cent surprising discoveries is self-sustained current oscil- 
lations (SSCOs) under a dc bias 3, 4]. A large number of 
experimental and theoretical studies have focused on dif- 
ferent aspects of these oscillations. It is understood that 
SSCOs are accompanied by the motion of boundaries of 
electric field domains (EFDs) Q . A model capable of de- 
scribing both the formation of stationary EFDs and SS- 
COs emerged after many tedious analyses and numerical 
calculations 0. Our understanding of SSCOs was greatly 
advanced through numerical investigations of this model. 
As we know, physical system with an intrinsic frequency 
may have various possible responses to an external driv- 
ing force. For a linear system, such as a simple pendu- 
lum, it will oscillate with the frequency of the driving 
force. For a nonlinear system, however, it can oscillate 
with an integer multiple of the driving frequency. For 
example, a laser light passing through a nonlinear opti- 
cal medium may lead to the second and third harmonic 
generations. It is also known that many other nonlinear 
systems can have frequency locking, in which a system 
might oscillate with an integer fraction of the driving 
frequency. It is therefore interesting to ask how the SS- 
COs observed in SLs response to a combined dc and ac 
bias. Numerical solutions on several SL models[^0,|^0 
show possible aperiodic oscillations, either quasiperiodic 
or chaotic. The experimental evidences of chaotic behav- 
ior were also reportedlH 113. While most early stud- 
ies focused on the chaotic behavior of tunneling current, 
there were also studies of the periodic response to an ac 
bias based on diople EFD model 12]. However, a gener- 
alized explanation with deep insights into the observed 
periodic response and other nonlinear responses in SLs is 
yet to emerge. 

In a recent study 0, Q|, we found that SSCOs in a 
SL correspond to the generation of limit cycles around 
an unstable steady-state solution. In the terminology 
of nonlinear physics0|, SSCOs are the manifestation of 
one-dimensional attractors-limit cycles. The power of the 



limit-cycle concept lies in its simplicity and universality. 
An important question one might ask is if the SSCOs 
indeed come from limit cycles, what will be the possible 
responses of the SSCOs to an extra ac bias? As it will 
be demonstrated in this letter, the EFDs model indepen- 
dent responses of SSCOs in SLs to an extra ac bias can be 
readily derived based on limit cycles. One will see that in 
periodic responses of SSCOs in SLs the frequency locking 
is not only a natural outcome of a limit cycle, but also 
very robust against the ac driving force. Furthermore, 
we can predict quantitatively the frequency of a periodic 
response for a given external ac bias. Thus, it is benefi- 
cial to use the limit cycle to understand SSCOs. First, 
by considering a limit cycle as a basic object, we argue 
that the limit cycle can have three possible responses 
to an extra ac bias: 1) a small deformation without a 
topological change; 2) a small deformation with a topo- 
logical change; 3) destruction of the limit cycle. The first 
two scenarios lead to the phenomenon of frequency lock- 
ing, and the last one gives rise to an aperiodic response. 
Second, a widely used drift velocity model is solved nu- 
merically to demonstrate this type of frequency locking 
and its robustness. Finally, we present our experimental 
results. An excellent agreement between the theory and 
the experiment is achieved. 

For the sequential electron tunneling in a SL, the phase 
space is made up by the bias on each potential barrier, 
as explained in References 13 and 14, since the state of 
a given SL is fully determined by these biases. On each 
point in the phase space, there is a unique vector which 
describes the system velocity in the space^^. This ve- 
locity is determined by the dynamics of the system. A 
point with zero velocity is called a fixed point. An unsta- 
ble fixed point, as denoted by the cross in Fig. 1, is such 
that a small deviation from the fixed point will drive the 
system away from the point. However, the system will 
stay around the fixed point because of the external bias 
constraint. In the case of a SSCO, this local repulsion and 
global attraction lead the system to move along a closed 
curve, a limit cycle, around the fixed point 0,0]. Using 
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a two-dimensional case as an example, it is schematically 
illustrated in Fig. 1 as a loop in solid line. Applying an 
extra ac bias with frequency ujac, the velocity field in the 
phase space changes through the dynamical equations of 
the system. If the bias is small, it can only perturb the 
velocity slightly, which, in turn, modifies the limit cycle. 

For the periodic response, a limit cycle can change in 
two distinct ways. One way is a small deformation of 
the limit cycle without a topological change, as shown 
by the curve A in Fig. 1. The length of the limit cycle 
can at most change a little. The time period for a sys- 
tem to move along the closed curve once does not change 
substantially, since its velocity field in the phase space 
is controlled by the system dynamics that is perturbed 
only slightly by the extra ac bias. In this case, the system 
oscillates with a frequency uj, whose value is not too far 
from its intrinsic frequency loq- There is another require- 
ment for the periodic response. Considering the system 
starting initially from a point on a limit cycle, it moves 
along the limit cycle and returns to the starting point 
after a time T, giving a frequency lo — 2t:/T. To have a 
periodic motion, the external ac bias should also return 
to its initial value. This means that T must be an inte- 
ger multiple of the ac-bias period 27r/wac- Thus, we have 
'^ac/'^ = P =integer. In fact, this is a general condition 
for the periodic motion of a dynamical system under an 
ac driving force^^. A natural conclusion of this argu- 
ment is that the limit cycle makes a small deformation 
like that of curve A in Fig. 1, when the ac frequency Woc 
is in the vicinity of puJo, an integer number of the system 
intrinsic frequency. 




FIG. 1: Schematic drawing of limit cycles around an unstable 
fixed point (cross) in a phase plane. The loop in solid line is 
the limit cycle in the absence of an ac bias. Curves A and 
B are the two possible deformation of a limit cycle under an 
ac bias. The system shall oscillate with a frequency close to 
its intrinsic one in case A while it oscillates with half of the 
intrinsic frequency in case B. 

The other way is a topological change of the limit cycle 
in spite of slight perturbation of the phase velocity field. 
This can occur when the system trajectory does not close 
itself after moving around the fixed point once. Instead, 



the trajectory returns to its starting point after q rounds. 
For example, the curve B in Fig. 1 shows a closed curve 
after two rounds. In this case, the system oscillates with a 
frequency around wo/2. It should be pointed out that this 
situation can occur only when the dimensionality of the 
phase space is larger than two. Combining the periodicity 
requirement discussed above, the second scenario occurs 
only when Wac/^o is in the vicinity of p/q integer), 
where p and q are integers. Under this type of ac bias, the 
limit cycle deforms itself in such a way that it becomes 
a closed curve after q turns in the phase space, giving 

UJac/l^ = P- 

Obviously, the periodic response must oscillate with 
an integer fraction of the ac frequency if the SSCOs are 
indeed due to the generation of limit cycles around an 
unstable fixed point. This type of frequency locking does 
not depend on a particular model. Since the limit cycle 
can deform itself under an ac bias, this periodic response 
with an integer fraction of the ac frequency is expected to 
be quite robust against tUac, meaning that ujac/i^ should 
remain unchanged with a small variation of ujac- In or- 
der to verify these results, we numerically solve a widely 
used discrete drift model under the combined dc and ac 
biases in the SSCOs regime. For a system consisting of 
N quantum wells under a bias U between the two end 
wells, the current flow is perpendicular to the SL layers. 
In the sequential tunneling, charge carriers are in local 
equilibrium within each well, so that a chemical potential 
can be defined locally. The chemical potential difference 
between two adjacent wells is called bias V on the bar- 
rier between the two wells. A current li passes through 
the i*^ barrier under a given bias Vi. This current may 
depend on other parameters, such as doping Nd. 

Following References 5 and 14, the dynamics of the 
system is governed by the discrete Poisson equations 

kiV,-V,^i) = n,- Nd, i^l,2,...N (1) 

and the current continuity equations 

dV 

J = k^+h, ?: = 0,1,2,...7V (2) 

where k depends on the SL structure and its dielectric 
constant, rii is the electric charge in the z*'' well. In 
Eq. the same doping in all wells is assumed, h is, 
in general, a function of Vi and n^. It can be shown|l7j 
that all SSSs are stable if li is a function of Vi only. On 
the other hand, a SSS may be unstable's'l if one chooses 
li — niv[Vi), where i; is a phenomenological drift velocity 
which is, for simplicity, assumed to be a function of Vi 
only. The constraint equation for Vi is 

N 

Y^V^U (3) 

To close the equations, a suitable boundary condition is 
needed. It is reasonable to assume a constant no, = 
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SNd , if the carrier density in the emitter is much larger 
than those in wells, and its change due to a tiny tunneling 
current is negligible. 

Previous studies [3, have shown that this model is 
capable of describing SSCOs with a negative differential 
drift velocity. One can obtain a SSCO when v{V) = 
0.0081/[(y/£; - 1)2 + 0.01] + 0.36/[(F/£; - 2.35)2 + 0.18] , 
iV = 30, J7 = 32.7E, No = O.OmkE, and S = 1.001 are 
usedQ. This V has two peaks atV ^ E and V = 2.35E. 
The region from V = E to V = 1.3E exhibits negative 
differential velocity. Thus, E can be used as a natural 
unit of bias, and l/v{E) as that of the time (the lat- 
tice constant is set to be 1). The intrinsic frequency loq 
is 0.14(w(i5)/l), indicating that the corresponding EFD 
boundary oscillates inside about 7 wells. Now we apply 
an extra ac bias Vac sin(waci) with Vac = 0.327E in ad- 
dition to the above dc bias. And we have solved numer- 
ically the above set of equations for different ujac- The 
current oscillation frequency can be obtained from the 
Fourier transformation of time evolution of the current. 
The results are plotted in the ujac/i^ vs. Uac/^o plane 
shown by the solid line in Fig. 2. It has a structure 
similar to a devil's staircase. The width of the staircase 
contains the information of robustness that the limit cy- 
cle can adjust itself. This robustness depends on both the 
amplitude and the frequency of ac bias. The width of the 
devil's staircase around Uac/t^o — ^/q decreases with the 
increase of q and the results are not displayed for g > 3. 
The dash line is the similar result for Vac — 0.654i?. The 
data are offset vertically for a better view. Clearly, the 
periodic response around Uac = LSwo disappears, lead- 
ing to a possible chaotic response. It shows that the limit 
cycle is destroyed, and the phase trajectory, somehow, 
cannot make a closed curve under this ac bias. 



response of SSCOs under an extra ac bias from a 
GaAs/AlAs SL sample. The GaAs/AlAs SL sample is 
grown by molecular beam epitaxy. It consists of 30 peri- 
ods of 14nm GaAs well and 4nm AlAs barrier and is sand- 
wiched between two n-l— GaAs layers. The central lOnm 
of each GaAs well is doped with Si (n = 2 x 10^ ''cm^'^). 
The sample is fabricated into 0.2 x 0.2mm2 mesas. The 
SSCOs response is recorded using an Agilent 54642A os- 
cilloscope. It has been found that SSCOs in a SL can 
be induced by changing the sample temperature^. In 
this measurement, the sample temperature is fixed at 
95K and the dc bias at 0.34V, which is located within 
the first plateau of the time-averaged I-V curve. The in- 
set of Fig. 3 shows the current oscillation trace without 
any ac bias. The SSCOs is clearly demonstrated with a 
frequency luq = 51.56KIIz. The oscillation periodicity is 
indicated by the corresponding Poincare map (or the first 
return map) as shown in Fig. 3. The Poincare map 
is derived from the current oscillation trace by sampling 
the current trace in a step of Tq (=27r/wo). Cares have 
been taken to minimize any possible artifact in deriving 
Poincare maps. 
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FIG. 2: ujac/io vs. ujac/ujo- The devil's staircase type be- 
havior with LUac/i^ = p, with integer p, shows that the system 
oscillates with an integer fraction of ac frequency, and this 
response is robust against ujac- Solid line is for Vac ~ 0.327i5 
and dashed line is for Vac = 0.654i5. Staircases with very 
narrow widths are not displayed. 

To verify the above results of frequency locking and 
its robustness experimentally, we have measured the 



FIG. 3: Poincare map: the points are {I(nTo), I(nTo 4- To)) 
for n = 1, 2, . . .. The accumulation of all points into a single 
dot in this map indicates that the current oscillates periodi- 
cally with a period of To . Inset: time-dependence of tunneling 
current. 

Figure 4 are Poincare maps obtained with an ap- 
plied extra ac bias. The ac bias amplitude Vac is set 
at 66mV and the driving frequencies ojac are indicated. 
The Poincare maps are obtained with sampling steps Tac 
(= 2'K/uJac)- As discussed early in this paper, the sys- 
tem exhibits the frequency locking when ujac is set in the 
vicinity of puj^ or LOa/p with p = integer. The measured 
response frequency uj is equal to uJac/p for uJac — P^o 
or LOac for LOac = ^o/p, respectively. Fig.4 (a) -(d) 
clearly demonstrate the occurrence of the frequency lock- 
ing when LOac — ojq, 2luo, luq/S, and ujq/2, respectively. 
The two dots in Fig. 4(b) result from uj = uJac/2 while 
our sampling step is Tac- 

In order to demonstrate the robustness of the fre- 
quency locking, we focus on the frequency locking in the 
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vicinity of ujac = 3a;o/2 and 2a;o- By varying the applied 
ac bias, the existence of the devil's staircase is shown in 
Fig. 5, in which the data for different ac bias amplitude 
Vac are offset vertically for clarity. Lines in Fig. 5 indicate 
the frequency locking range around Wac/wo = 1.5 and 2 
for a given Vac- Clearly, the devils staircase widths, i.e. 
the robustness of frequency locking, are strongly depen- 
dent on Vac and Uac- The locking range for u>ac/^o = 2 
is much larger than that for LOac/^jJo ~ 1-5. When Vac 
is small (=16mV) the periodic response cannot be found 
around uiacl'^^^ = 1-5. These results are in good agree- 
ment with the theoretical expectations. 

The above demonstrated quantitative agreements 
among the theory, the model calculation, and the ex- 
periment clearly indicate that SSCOs in SLs are indeed 
originated from the generation of limit cycles. The un- 
usual frequency locking into a particular set of an integer 
fraction of ac frequency is the direct manifestation of de- 
formation of limit cycles under ac bias. It is worth to 
emphasize that based on limit cycles frequency locking 
does not depend on the particular structure of the EFD. 
Like many other nonlinear dynamical systems, the fre- 
quency locking of SSCOs can be understood within the 
framework of the general concepts and principles of non- 
linear physics. 

In summary, the limit cycle can deform itself in such 
a way that it makes q turns in the phase space around 
an unstable fixed point when the ac frequency LOac is in 
the vicinity of ^wq with integers p and q. Thus, a system 
may oscillate with Uac/p, an integer fraction of LOac, or 
about Wo /q- Both of this type of frequency locking and its 
robustness are verified by the numerical model calcula- 
tions and real experiments. In conclusion, the limit cycle 
picture of SSCOs gives a deep insight into the nonlinear 
properties of SLs. 
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FIG. 4: Poincare maps with sampling step Tac for cjac ~ i^o 
(a); Wac = 2uJo (b); Wac = tJo/3 (c); ujac = t^o/2 (d). The num- 
ber of dots in the maps multiplying Tac are the corresponding 
response periods. Thus, the frequencies for (a), (c) and (d) 
are all equal to Uac while the frequency for (b) is Uacfi- 



FIG. 5: Experimentally measured periodic response ranges 
around ujac/i^o ~ 1-5 and 2, respectively, illustrated in ujac/i^ 
vs. uiac/i-oo plots for difi'erent amplitudes of ac bias indicated. 
The data are offset for clarity. 



